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ABSTRACT 

We have obtained optical and near-infrared images of the field of the accreting mil- 
lisecond X-ray pulsar XTE J1751-305. There are no stars in the 0.7" error circle (0.7" 
is the overall uncertainty arising from tying the optical and X-ray images and from 
the intrinsic uncertainty in the Chandra X-ray astrometric solution) . We derive limit- 
ing magnitudes for the counterpart of R>23.1, I>21.6, Z>20.6, J>19.6, K>19.2. We 
compare these upper limits with the magnitudes one would expect for simple models 
for the possible donor stars and the accretion disk subject to the reddening observed 
in X-rays for XTE J1751-305 and when put at the distance of the Galactic Centre 
(8.5 kpc). We conclude that our non-detection does not constrain any of the models 
for the accretion disk or possible donor stars. Deep, near-infrared images obtained 
during quiescence will, however, constrain possible models for the donor stars in this 
ultra-compact system. 
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1 INTRODUCTION 

Low-mass X-ray binaries (LMXBs) are typically old 
(> 10 8 yr) binary systems in which a low-mass companion 
star (^IMq) transfers matter to a neutron star or a 
black hole. The neutron star LMXBs are thought to be 
among the predecessors of the millisecond radio pulsars. 
Due to accretion of matter and decay of the magnetic 
field during the LMXB phase, the neutron star spins-up 
to millisecond periods (see Bhattacharya 1995 for a re- 
view). However, until March 2002 only six LMXBs were 
known to show pulsations (see IJonker fc van der Kli s 200T1 

* email : peterj@ast.cam.ac.uk; based on observations collected 
at the European Southern Observatory, Chile (ESO No 269. D- 
5013) 



for a short overview), and only one of them, a tran- 
sient system, was shown to have a millisecond period 
(SAX J18 08.4-3658; |Wfjn"ands fc van der Klis 1998] 
and |Chakrabarty <£TTyio7glm"l*9*9*8f ! In 

2002 |Markwardt fc Swank (2002) | (see also 
IMarkwardt et al. 20021 and |Galloway et al. (2002b71 
(see also Gallowa y et al. 2002at discovered millisecond 
pulsations in two other transient sources, XTE J1751-305 
and XTE J0929-314, in outburst. 

There are two main reasons to search for the compan- 
ion stars of accretion powered millisecond X-ray pulsars. 
Accretion powered millisecond X-ray pulsars have most 
likely accreted a substantial amount of matter as they 
were spun-up by accretion to a period of milliseconds. A 
spectroscopical determination of the radial velocity curve 
of the companion star yields a lower limit to the mass 
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of the neutron star; measuring a mass of considerably 
more than 1.4 Mq for even one neutron star would imply 
the firm rejection of many proposed equations of state 
(see the discussion by |van Paradijs fc McClintock 1995| 
for an overview of the mass determinations of millisecond 
radio pulsars in neutron star-white dwarf binaries see 
|Thorsett fc Ch akrabarty 19990 . In order to spectroscopi- 
cally measure the radial velocity of the companion star, 
clearly a detection of the object is first required. Secondly, 
at present there are three known ways to produce compact 
X-ray binaries (Savonije, de Kool 



of XTE J1751-305 on the CCD in our image obtained with 
the Danish telescope, rendering it impossible to determine 
whether a source is present or not at the Chandra X-ray 
location. Hence, we will not consider this observation any 
further. 

We derived an astrometric solution for the optical 
I band image obtained with the NTT of the field of 
XTE J1751-305 using the positions of 4 nearby unsaturated 
stars which appear in the USNO-A1.0 catalogue. The rms 
of the fit was 0.015". The typical astrometric error of stars 
in the USNO-A1.0 catalogue is 0.25" (68 per cent confi- 



^^^nderiF£ewel^86l 

I Tutukov fc Fedorova 1989l | Y u ngelson, JNelemans fc van den Heuv^»^| ^ c uncertainty in the X-ra7position is dominated 
ITutukov et al. 1987l|Pod siadlowski, Rappa port fc Pfahl 2002| 



see Section 3 for a more elaborate account of these models) . 
Detection of the companion star will allow us to differentiate 
between these different formation scenarios. 

In this Paper we present the results of our search for the 
optical and near-infrared counterpart of the ultra-compact 
accreting millisecond X-ray pulsar XTE J1751-305 
(P orb = 42 minutes; IMarkwardt et al. 20*021 . Early 
reports on the data presented in this Paper were 
published by |Jonker et al. (2002)] fWang et al. (2002)| 
|Wang fc Chakrabarty (2002)| and |Kong et al. (2002)| 



2 OBSERVATIONS AND ANALYSIS 

We observed the region around the accurate Chandra X- 
ray position of the source using the 6.5 m Magellan, the 
3.8 m UKIRT, the 3.58 m NTT, and the 1.54 m Danish 
telescopes. A spectrum of the candidate optical counterpart 
presented by |Kong et al. (2002)| was obtained with ESO's 
3.6 m telescope at La Silla using Grism number 12 which 
covers 601-1032 nm. A log of the observations can be found 
in Table In case of the near-infrared data the total ex- 
posure time is given and the dither pattern, i.e. in case of 
the UKIRT J band images we took 5 exposures of 60 sec- 
onds each, in between each observation the telescope was 
dithered; this pattern was repeated twice. In the case of the 
Magellan near-infrared observations we obtained J and Ks 
images of 60 seconds integration time each. We also obtained 
11 J band images of 20 seconds integration each. 

The data reduction was performed in iraf 1 . The near- 
infrared images (J and K band) were sky subtracted (using 
the sky determined from the dithered images), flatfielded, 
aligned, and combined to form one image per band per night. 
The Magellan near-infrared data have been linearised before 
the reduction. The optical images (R, I, and Z band) were 
reduced in the standard fashion. Aperture and point spread 
function fitting photometry were done using the packages 
APPHOT and daophot in iraf. The spectrum obtained with 
ESO's 3.6 m telescope was extracted using the task APALL 
in the package SPECRED in iraf. All nights were photomet- 
ric except for the nights of April 13 and 14, 2002 when the 
humidity was very high at Magellan and April 14, 2002 at 
UKIRT when thin cirrus was present. The J band observa- 
tions during that night also suffered from a bright sky as 
they were obtained partially during twilight. Charge over- 
flow from a neighbouring bright star covered the location 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories 



by the Chandra aspect solution (0.6" ; 90 per cent confidence; 
IMarkwardt et al. 200211 . Hence, the overall astrometric un- 
certainty in the position of the errorcircle is 0.7" (90 per cent 
confidence). Next, we tied the astrometric solution of the I 
band image to the optical R and Z band and to the near- 
infrared J and K(s) band images by assigning the known 
position of several stars (from the I band) to the R, Z, J(s), 
or K band image; this did not increase the error in the as- 
trometry significantly. 

In Figure {top panel) we show the R band image 
from our NTT observations. Separately, in Figure (bot- 
tom panel) we show the near-infrared UKIRT (K band) 
image. From the astrometry it is clear that each of the 
two candidate counterparts (both indicated with two ar- 
rows in the top panel) is just outside the 90 per cent confi- 
dence error circle ( |Wang et al. 2002 | initially proposed star 1 
as the counterpart but see also |Wang fc Chakrabarty 20021 
Kong ~ al. 2002l proposed star 2 as the counterpart; see Fig- 
ure 0. Hence, it is unlikely but not ruled out that the coun- 
terpart was detected. 

We determined 3 o upper limits on the presence of a 
star in the error circle of each of the images obtained under 
photometric conditions. To do this we added a simulated 
artificial star (created using the point spread function of the 
stars in the image) at the position of the error circle. The 
magnitude of this extra star was varied and measured using 
the standard photometric tasks. We define the 3 o limit- 
ing magnitude as the magnitude at which the error in the 
magnitude of such an artificial star is ~0.3 magnitudes (an 
0.3 magnitude error is equivalent to a ~30 per cent error 
on the flux measurement, i.e. a three sigma detection; see 
also Hullc man et al. 200011 . For all the images we need ~3-5 
trials to cover the error circle (the number of trials is dif- 
ferent for the different images; it varies as a function of the 
seeing). We did not take the number of trials into account 
when calculating the upper limits (the error introduced by 
this is ~0.1 magnitude). The error on the determined lim- 
iting magnitude is also ~0.1 magnitude. Together with the 
uncertainty in the photometric zero-point (at most 0.1 mag- 
nitude) this yields an uncertainty of ~0.2 magnitude in the 
limiting magnitudes. 

For the NTT Z band image we converted the (Sloan) z' 
magnitude of the observed standard star (PG0918+029D) 
to the UKIRT UFTI Z band using the transformation given 
on the UKIRT web page. However, since the NTT and the 
UKIRT Z band differ, and there is no cross-calibration of 
these two Z bands, in doing so we introduced an uncer- 
tainty in the Z band magnitudes which can be more than 
0.2 magnitudes. Therefore, the Z band upper limit is given 
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Table 1. Log of the observations. MJD and start time refer to the time at the start of the first observation. 



Telescope 


Instrument 


Observation date 


MJD 


Filters 


Exposure 


Airmass 


Seeing 






and start time (UT) 


(UTC) 


Bessel (except J, K(s)) 


Time (s) 




(arcseconds) 


Danish 


DFOSC 


08-04-2002 08:51 


52372.36870 


R 


900 


1.01 


~1.2 


Magellan 


Classic-CAM 


09-04-2002 08:02 


52373.33472 


Ks 


60 


1.05 


~0.65 


Magellan 


Classic-CAM 


09-04-2002 07:32 


52373.31388 


J 


280 (11x20+60) 


1.05 


~0.65 


Magellan 


MagIC 


13-04-2002 07:43 


52377.32153 


1 


240 (3x) 


1.05 


~0.6 


Magellan 


MagIC 


14-04-2002 08:30 


52378.35416 


I 


30 (lx), 240 (18x) 


1.00-1.03 


~0.7 


UKIRT 


UFTI 


14-04-2002 14:41 


52378.61201 


.1 


600 (2x5x60) 


1.6 


~0.8 


UKIRT 


UFTI 


14-04-2002 14:55 


52378.62161 


K 


600 (2x5x60) 


1.6 


~0.65 


UKIRT 


UFTI 


18-04-2002 14:47 


52382.61623 


.1 


600 (2x5x60) 


1.6 


~0.65 


UKIRT 


UFTI 


18-04-2002 15:06 


52382.62925 


K 


600 (2x5x60) 


1.6 


~0.5 


NTT 


SUSI2 


18-04-2002 08:33 


52382.35682 


R 


600 


1.00 


~0.8 


NTT 


SUSI2 


18-04-2002 08:23 


52382.34950 


I 


600 


1.00 


-0.9 


NTT 


SUSI2 


18-04-2002 08:12 


52382.34185 


Z 


600 


1.01 


-0.8 


3.6M (spec) 


EFOSC 


01-05-2002 07:19 


52395.30470 


Gr#12 0.7" slit 


1500 


1.01 


~0.9 



Table 2. Upper limits (3 <r) on the presence of a star at the position of the Chandra error circle. 



Telescope Observation date Limiting magnitude 



Magellan April, 9, 2002 J >19.6 

Magellan April, 9, 2002 K >18.2 

UKIRT April 18, 2002 J >19.5 

UKIRT April 18, 2002 K >19.2 

NTT April 18, 2002 R>23.1 

NTT April 18, 2002 I>21.6 

NTT April 18, 2002 Z>20.6 



as reference only; it should be considered approximate. The 
upper limits are given in Table [5] 

We searched for variability in the I band mag- 
nitude for the candidate counterpart proposed by 
|Wang fc Chakrabarty (2002) | (star 1) by comparing the 
magnitudes in the Magellan and the NTT observations. 
The I band magnitude of star 1 is consistent with be- 
ing the same during the observations. Unfortunately, the I 
band magnitudes of the candidate counterpart proposed by 
|Kong et al. (2002)| (star 2) could not be determined for the 
NTT image due to the presence of the nearby star. However, 
differential photometry showed that star 2 was not variable 
in the Magellan Sloan i' band images taken a day apart 
HKong et al. 2002| . 

The slit we used to obtain a spectrum using ESO's 3.6 m 
telescope at La Silla had a width of 0.7". The slit orientation 
was such that both star 1 and 2 were in the slit as well as 
most of the position marked by the error circle. However, 
due to the seeing of ~0.9" during these observations the 
spectrum of star 2 may have been contaminated by light of 
the bright nearby star. In order to investigate this we con- 
verted the magnitudes of this bright nearby star (R=18.3, 
1=16.2) and those of star 2 (R=22.6, 1=19.2; we note that 
the uncertainty on these magnitudes is large due to the non- 
photometric conditions at the time of the Magellan I band 
observations; for comparison the best estimate of the mag- 
nitude of star 1 was 20.2) to flux densities and we compared 
the amplitude of the spectral energy distribution according 
to these R and I band flux densities with the amplitude of 
the flux calibrated spectrum. The fact that the flux density 



of the spectrum is higher than that of star 2 shows that 
the spectrum will have been contaminated with light from 
the bright nearby star. We note that in converting the mag- 
nitudes to fluxes we assumed that star 2 was not variable. 
We show the spectrum in Figure [5] We also label the po- 
sition of Ha, although it is unclear whether the donor star 
of XTE J1751-305 contains hydrogen or not. The spectrum 
is featureless except for the atmospheric absorption feature 
near 7613 Angstrom and the band at ~9300 Angstrom. The 
spectrum of star 1 was not detected, rendering further sup- 
port to the conclusion that the detected light from the po- 
sition of star 2 was dominated by the bright nearby star. 



3 DISCUSSION 

We have obtained optical and near-infrared images of the 
field of the accreting millisecond X-ray pulsar XTE J1751- 
305. Two possible counterparts have been investigated but 
since they both fall outside the 90 per cent confidence Chan- 
dra error circle we conclude that neither the optical nor 
the near-infrared counterpart was detected. The star clos- 
est to the error circle was suggested to be the counterpart 
by |Kong et al. (2002)| (star 2 in Figure[TJl. We placed upper 
limits on the presence of a star in the error circle in the R, 
I, Z, J, and K band. 

To evaluate the constraints these upper limits can place 
on the properties of the system, we consider the possi- 
ble evolutionary states for XTE J1751-305. To arrive at 
an approximately forty minute orbital period X-ray bi- 
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Figure 1. Top panel: The R band image (~ 27" X 14"; 10 minutes integration; North is up, East to the left) of the region of XTE J1751- 
305 obtained with the 3.58 m NTT. The X— ray error circle for the location of the source is overplotted (0.7"; 90 per cent confidence). 
The double set of arrows indicate the location of the possible counterparts proposed by Wang & Chakrabarty (2002; star 1) and Kong 
et al. (2002; star 2). Bottom panel: The UKIRT K band image obtained on April 18, 2002 (~ 27" X 14"; 10 minutes integration; North 
is up, East is to the left) of the region of XTE J1751— 305. The Chandra error circle is overplotted. The vertical black strips in the centre 
of the image are artifacts introduced when combining the dithered images. 



nary there are essentially three evolutionary paths. The 
first starts from a detached white dwarf - neutron star 
binary, which is brought into contact by angular momen- 
tum loss due to gravitational wave radiation. For a discus- 
sion of this scenario for ultra-compact X-ray binaries see 
|Yungelson, Nelemans fc van den Heuvel (2002)| The mass 
of the donor in this case would be ~0.02 Mq and the mass 
transfer rate ~10" n M yr _1 . Recently, |Bildsten (2002"jl dis- 



cussed XTE J1751-305 and the similar accreting millisecond 
X-ray pulsar XTE J0929-314 and calculated models for hot 
white dwarf donors, rather than the previously used zero- 
temperature models. He finds a luminosity of the order of 
1CP 3 L.0 and an effective temperature of ~6000 K for the 
donors. 

It is also possible to form a forty minute binary 
from a helium star that transfers matter to a neu- 
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6000 7000 8000 9000 1 O 4 " 

Wavelength (Angstrom) 

Figure 2. Spectrum of the Southernmost source identified with the marks in Figure Q (star 2; to the East of the bright neighbour), 
obtained with ESO's 3.6 m telescope at La Silla on May 1, 2002. The line connects the magnitudes of this star in the R and I band 
converted to flux densities at the central wavelength of the R and I band filter. 



tron star (e.g. |Savonije, de Kool fe van den Heuvel 1986| 
ITutukov fc Fedorova 198911 . Such' a scenario goes through 
the forty minute period range twice. The first time, the sys- 
tem goes through the orbital period of 40 minutes while the 
orbital period decreases, then the donor is a luminous, hot 
helium star of about 0.6 Mq. The mass transfer rate at that 
stage is of the order of 10 -8 M yr -1 . After having reached 
a period minimum of about 10 minutes, the system returns 
to longer periods with a dim, semi-degenerate helium star 
donor, not unlike the low-mass white dwarfs in the first sce- 
nario. 

The third evolutionary scenario producing ultra- 
compact binaries involves a main sequence star 
close to core hydrogen exhaustion that starts mass- 
transfer to a neutron star {Tu tukov et al. 19871 
|Podsiadlowski, Rappaport & Pfahl 2002). Such a sys- 
tem evolves to shorter orbital periods than the standard 
period-minimum for hydrogen rich stars, since more 
and more helium enriched layers are reached when peel- 
ing off the star. These systems go through a period 
minimum, again of the order of 10 minutes, and form 
white dwarf like, low-mass, low-luminosity donor stars. 
|Podsiadlowski, Rappaport fc Pfahl (2002) | discuss possible 
models for XB 1832-330 which has an orbital period close 
to forty minutes as well, so these models should also be 



applicable to XTE J1751-305. The donors in these models 
have masses and mass transfer rates between 0.026 and 
0.14 M and 4.5 and 22xl0" n Mq yr -1 , respectively. 

The X-ray observations already place some con- 
straints on the system. The mass-function for the com- 
panion (1.3xl0~ 6 M ; IMarkwardt et al. 2 002 1 suggests 
that the companion is a low-mass star. If the companion 
were to be a 0.6 Mq helium star, the inclination should 
be less than 1.6°. The two pre period-minimum mod- 
els of |Podsiadlowski, Rappaport fc Pfahl (2002)| with donor 
masses of 0.14 and 0.094 M would imply inclinations less 
than 6° and 9° , respectively. |Markwardt et al. (2002)] infer 
a mass transfer rate of 2. 1 x 10~ Mq yr~ , assuming a neu- 
tron star radius of 10 km, a distance of 10 kpc and a recur- 
rence time of 3.8 yr. Even though there is quite some un- 
certainly in the distance and recurrence time, the high mass 
transfer rate as expected from a helium star donor seems 
to be unlikely, even more so, because at such high rates the 
system would be expected to be a persistent X-ray source 
HTsugawa fc Osaki 19971 - 

Next, we estimate the expected absolute magnitudes of 
the different donor star and accretion disk models in the 
optical and near-infrared. Since the source was still in out- 
burst when we obtained our observations it is likely that the 
accretion disk dominates in the optical and near-infrared. 
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Figure 3. Top panel: Absolute magnitudes for the theoretical 
models for the donor stars in 40 minute periods ultra-compact 
X-ray binaries. Solid line: the hot white dwarf model from 
|Bildsten (2002) | Dotted lines: the pre period-minimum models 
for XB 1832—330 of (left two triangles in the top panel of fig- 
ure 16 of |Podsiadlowski , Rappa port &: Pfahl 2002*} . Dashed lines: 
estimates for the outburst disc spectrum, based on an absolute 
visual magnitude of 4.5 and temperatures of 10000 and 30000 
K. The top dashed line also gives the absolute magnitudes of 
the donor in case it would be a luminous, hot helium star (from 
Tutukov & Fcdorova 19891. The passbands we used for the filters 
are indicated (we assumed perfect step function filters). Bottom 
panel: The estimated magnitudes for the same models as in Fig- 
ure [3] assuming the system is close to the Galactic centre and 
taking the reddening into account (Ay = 5.6, based on the esti- 
mated Njj ~ 1 X 10 22 cm -2 ). The filled triangles represent the 
upper limits on the counterpart of XTE J 175 1-305. 



absolute visual magnitude. Since these ultra-compact sys- 
tems have to have small accretion disks, they are expected 
to be hot (see, |van Paradijs fc McClintock 1995| l. We esti- 
mated the absolute magnitudes for disks modelling them as 
simple blackbodies of 10000 and 30000 K. The results are 
shown in Figure [2] (dashed lines in the top panel). 

The absolute magnitudes of the possible donor 
stars in a forty minute binary, again assuming sim- 
ple blackbody spectral energy distributions, are calcu- 
lated for the hot white dwarf scenario proposed by 
|Bildsten (2002)|and for the two pre period-minimum mod- 
els of |Podsiadlowski, Rappaport fc Pfahl (2002)1 We did not 
consider their post-period minimum model because accord- 
ing to |Bildsten (2002)1 the companions in X-ray binaries will 
never cool down so much. Also, such a donor star would be 
very similar to the hot white dwarf donor model. For the hot 
helium star donor, we used a model having a luminosity of 
100 L Q and T = 60000 K (see lTutukov fc Fedorova 198911 . 
Interestingly, the absolute magnitudes of this model fall on 
top of the 30000 K disc model. The absolute magnitudes of 
all four models are shown in Figure [3] (top panel). 

To compare the theoretical models with our upper lim- 
its, we assumed XTE J1751-305 is near the Galactic Cen- 
tre, with a distance modulus of 14.65 magnitudes. Fur- 
thermore, we inferred an absorption (Ay — 5.6), based 
on the measured Nh = 10 cm - UMiller et al. 20021 and 
calculated the absorption in the other bands according to 
the relations found by |Rieke fc Lebofsky (1985")] The re- 
sults are shown in Figure [3] (bottom panel) where the sin- 
gle symbols denote the upper limits. Unfortunately, our 
upper limits do not constrain any of the models. Fur- 
ther deep near-infrared imaging could, however, start to 
rule out several possible models for the counterpart of the 
accreting millisecond X-ray pulsar XTE J1751-305. Ac- 
cording to the models presented above, XTE J0929-314, 
having an orbital period of ~40 minutes, a reddening in 
the V band of ~0.65 magnitudes (Dic key fc Lockman 1990} 
can't have a distance much larger than the lower limit 
of 6.5 kpc derived by |Galloway et al. (2002a) | Otherwise 
it would not have been detected in outburst at V=18.8 
iGreenhill, Giles fc Hill 2002} . 
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